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BINARITY IN COOL ASYMPTOTIC GIANT BRANCH 
STARS: A GALEX SEARCH FOR ULTRAVIOLET EXCESSES 

R. Sahai^, K. Findeisen^, A. Gil de Paz"^, C Sanchez Contreras^ 

sahai® jpl . nasa . gov 

ABSTRACT 



The search for binarity in AGB stars is of critical importance for our un- 
^ I derstanding of how planetary nebulae acquire the dazzling variety of aspherical 

shapes which characterises this class. However, detecting binary companions in 
' r-|'' such stars has been severely hampered due to their extreme luminosities and pul- 

sations. We have carried out a small imaging survey of AGB stars in ultraviolet 



P I light (using GALEX) where these cool objects are very faint, in order to search 



H 

c/3 I for hotter companions. We report the discovery of significant far-ultraviolet ex- 

cesses towards nine of these stars. The far-ultraviolet excess most likely results 
either directly from the presence of a hot binary companion, or indirectly from 

P' \ a hot accretion disk around the companion. 



Subject headings: binaries: general, planetary nebulae: general, stars: AGB and 
post-AGB, stars: mass-loss, circumstellar matter 



1. Introduction 



X 

Vh ' There are many observational indications which lead us to believe that binarity, believed 

to be very common amongst pre- main-sequence (e.g., Bodenheimer et al. 2000) and main- 
sequence stars (Duquennoy & Mayor 1991), strongly influences the history and geometry of 
mass loss during the late stages of stellar evolution. The evolutionary transition from the 
AGB to the post-AGB phase is accompanied by significant changes in the morphology of 



iJet Propulsion Laboratory, MS 183-900, Caltech, Pasadena, CA 91109 

^Cornell University, Ithaca, NY 14853 

■^Dpto. de Astrofisica, Universidad Complutcnse de Madrid, Madrid 28040, Spain 

^Dpto. de Astrofisica Molecular e Infraroja, Instituto de Estructura de la Materia-CSIC, Serrano 121, 
28006 Madrid, Spain 



-2- 



these objects - the roughly round circumstellar mass-loss envelopes (CSEs) of AGB stars 
evolve into post- AGB nebulae with a dazzling variety of shapes and intriguing symmetries 
(e.g. Schwarz et al. 1992, Sahai & Trauger 1998). Critical reviews (Soker 1998) of the 
properties of bipolar PNe (e.g. Corradi & Schwarz 1995) lead to the conclusion that binary 
models can explain all these properties, whereas single-star models (e.g. Garcia-Segura 1997) 
have many difficulties. 

However, in spite of dedicated efforts by many researchers to search for binarity in 
evolved stars, direct observational evidence for binarity has been hard to come by. AGB 
stars are very luminous (~ few x 10^-10^ L©) and surrounded by dusty envelopes, making 
it very difficult to directly detect nearby stellar companions which are generally likely to 
be significantly less luminous main-sequence stars or white dwarfs. Indirect techniques such 
as radial-velocity measurements (e.g., van Winckel et al. 1999, Sorensen & Pollaco 2003, de 
Marco et al. 2004) or photometric variability measurements (Bond 2000) have been used for 
the central stars of PNs and post-AGB objects, with some success. But these techniques 
cannot be easily applied to AGB stars, because the latter show strong variability intrinsic to 
their pulsating atmospheres, which potentially masks the corresponding variability due to a 
companion. Extensive observations of the central stars of planetary nebulae have resulted in 
detections a sum total of <20 binaries (Bond 2000, CiarduUo et al. 1999), implying a 10-15% 
fraction of detectable close binaries among randomly selected PNe. Bond (2000) concludes 
that it is likely that the known short-period binaries in PNe are only the tip of an iceberg 
of a substantial population of longer-period binaries. 

Deep ultraviolet observations hold the promise of being able to discover substantial 
numbers of binary companions in AGB stars, since most mass-losing AGB stars are relatively 
cool objects (spectral types ~M6 or later). The companions are likely to be main-sequence 
stars because of the steep dependence of evolutionary rates on stellar mass (e.g. Soker & 
Rappaport 2000). Thus for a secondary-to-primary mass ratio, q = M2/M1, around unity, 
any stellar companion has a good probability of being hotter than the primary. However, 
it is difficult to estimate with confidence the number of such systems where the secondary 
is on the main-sequence and hotter than the primary, as a fraction of the total number of 
primordial binaries, since the mass-ratio probability distribution, f{q), where q = M2/M1, 
is not well-known. A promising approach is to carry out population synthesis studies (which 
are still in their infancy) such as those of Soker & Rappaport (2000), who adopt, for their 
modelling, /(g) oc q^^^; note that this function is not strongly peaked towards q = I- 

Since observed and model spectra of cool AGB stars show that their fluxes die rapidly 
at wavelengths shortwards of about 2800A, signiflcantly favorable secondary-to-primary flux 
contrast ratios (>10) for companion detection may be reached in the GALEX FUV (1344- 
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1786A) and NUV (1771-2831A) bands, for companions of spectral type hotter than about 
GO (Te/j=6000 K). In this paper, we report on a subsample of objects from our Cycle 1 pilot 
program, which were detected in both the FUV and NUV bands, and the implications of 
these detections for binarity. A comprehensive study covering the full results of our survey 
will be presented in a forthcoming paper. 



2. Observations &; Results 

We selected a sample of 25 AGB late-M (i.e., M5 or later) stars (which passed the 
GALEX mission "bright-star" and "high-background" tests) largely based on their inclusion 
in the Hipparcos astrometric catalog, with a "multiplicity" flag in the header field H59 of 
the main catalog, indicating that a single-star astrometric solution was not adequate). Thus 
for 20/25 objects, the selection criteria of our "pilot" program were intentionally biased 
towards optimising the a priori probability of finding companions, in order to test the validity 
of our technique. Three objects did not have a "multiplicity" flag; and two are not in 
the Hipparcos catalog, but were selected from published lists of AGB stars with molecular 
envelopes detected in CO emission. Although most of our objects had positive entries for 
the annual parallax, the errors were usually large, and only for four objects were the parallax 
measurements significant (i.e., greater than 3a). The requirement that our objects be M5 
(or later) or cool N-type carbon stars was included in order to minimise the ratio of their 
UV fluxes to that of hotter companions (if present). 

From our original list of 25 objects, 21 objects have been observed - 20 as part of our 
GI program (GIl-23; PI: R.Sahai), and 1 as part of other programs - in both the NUV 
(1771-2831A) and FUV (1344-1786A) bands. Nine of these were detected in both the FUV 
and NUV bands with high S/N (> 8a) (Tabled . Amongst these, the NUV image of AFPeg 
has an elliptical shape, and an intensity cut along the major axis (PA~ 75°) shows two peaks 
separated by about 6"; the stronger one corresponds to AF Peg's location; in the FUV image, 
AFPeg is weaker. Because the separation is comparable to the PSF, the two sources cannot 
be deconvolved reliably; hence the measured fluxes of AF Peg are very uncertain. We have 
used the pipeline-generated catalogs included with the imaging datasets for extracting the 
photometry for the remaining 8 sources. Stars not detected in the FUV, will be discussed 
together with their NUV properties, in a forthcoming paper. 



^these were the so-called problem stars for which a double star solution with the classical parameters of 
separation and position angle could not be found: and were flagged as G, O, V, or X in field H59 of the main 
catalog 
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The 8 FUV sources include 4 oxygen-rich stars (RWBoo, AACam, VEri & RUMa) 
and 4 carbon-rich stars (TDra, TWHor, VHya, & VYUMa). For those objects where more 
than one exposure was available, taken at different epochs, we list both the individual and 
average ffuxes. The typical uncertainty in the measured ffuxes is dominated by systematic 
uncertainties in the GALEX pipeline photometric calibration of about 10-15% (Morrissey et 
al. 2005). 

We now consider whether our FUV and/or NUV detections could result from the pres- 
ence of a small filter red leak, which could produce spurious detection of a UV signal for the 
extremely red stars observed in this study. Since the GALEX detectors are photon counting 
MAMA detectors, there is supposedly no "red leak", since only UV photons can trigger 
the photoelectrons (Rich 2005); the photocathode on the FUV detector is non-responsive 
above ~1800A, and the NUV response is suppressed below measurable levels by multilayer 
coatings on the optics. According to the GALEX helpdesk, there is no measurable red leak 
in either GALEX band. 

Even though no red leak response has been measured for the GALEX filters, we have 
ensured that even if such a response is present at a low level, our modelling is not affected 
by the latter because the upper limits which we can set on the red leak from our data are 
quite low. We have done this by comparing, for our survey objects, the ratios of the FUV 
and NUV fluxes to the V-band fluxes and assuming that the lowest of these ratios are due 
to a red leak. From this analysis, we find values of 2.5 x 10~^ and 4 x 10~^ for the maximum 
possible "red leak" flux in the GALEX FUV and NUV bands, as a fraction of the V-band 
flux - these ratios are too low to affect our models, or the detection statistics we report in 
this paper. 



3. Ultraviolet Excesses 

We now investigate the origin of the ultraviolet fluxes in the objects we have detected 
in the FUV band. We have fitted the spectral-energy-distributions (SEDs) from 0.1 to 
2/im(Fig{T]) of the four oxygen-rich stars with reliable FUV fluxes (i.e. RWBoo, AACam, 
VEri and RUMa), using stellar atmosphere models of AGE stars (Fluks et al. 1994), 
corresponding to the spectral type of each star as given in the General Catalogue of Variable 
Stars (GCVS: Samus et al. 2004j). A visual extinction. Ay, to account for the extinction 
by circumstellar dust due to the dusty mass-loss envelope of the primary, is also determined 
from our fits. Archival photometry at wavelengths redwards of the GALEX NUV band was 



^available online through the Vizier astronomical database 
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taken from the Hubble Guide Star Catalog (GSC 2.2), the US Naval Observatory USNO-Bl.O 
Catalog, and the Two Micron All-Sky Survey. 

We find that the observed FUV (NUV) fiuxes are a factor > 10^ (>5) larger than 
expected for the photospheric emission of the primary, accounting for the finite filter band- 
width, the filter responsqfl, and the steeply sloping spectrum in the UV. Hence, even though 
the photometric variability of the primary stars makes our fit somewhat uncertain, it cer- 
tainly cannot account for the FUV excesses, because they are very large, and there is no 
systematic relationship between the light-cycle phases of the various photometric data-points 
used to fit the primary model which could conspire to produce such an excess in each of our 
three sources. 

The detailed SED fitting described could not be carried out for the 4 carbon-rich stars 
in Tabled] because, for these objects, model atmospheres for A < 2300A are not available 
(D. Luttermoser, priv.comm.). We therefore used black-body spectra using Te// values from 
Bergeat et al. (2001), and scaled these to fit the NIR and optical photometry of each object. 
We find that the model FUV fiux of the primary AGB star, even for the object with the 
smallest FUV excess (VYUMa), is lower than the observed value by a factor ~30. VHya 
stands out amongst the 9 FUV-detected stars as having the highest FUV-to-NUV fiux ratio 
(> 1), but the coolest photosphere (Tg// = 2160 K) for the primary. We discuss our detection 
of the FUV/NUV excess in this object in more detail later (§ll]). 

We note that Wood & Karovska (2004) conclude, based on lUE spectroscopic data of 
several Mira AGB stars at multiple phases, that this class of objects do not produce any 
detectable emission below 2000A. We have examined the lUE database for the sources in 
our survey, and find that only 3 objects in our survey sample were observed - R LMi, V 
Hya and TW Hor. R LMi and TW Hor were observed with both the long-wavelength (LWR 
and LWP) and short- wavelength (SWP) instruments (Boggess et al. 1978a,b), whereas V 
Hya was observed only with the long- wavelength instrument. For both R Lmi and V Hya, 
the signal-to-noise ratio over the observed bandpasses (1910-3300 A for LWR/LWP and 
1150-1975 A for the SWP) is close to zero. For TW Hor, the situation is similar, except 
at wavelengths longer than ~2500A, where significant fiux is detected; the steady rise in 
the spectrum towards the red end of the bandpass indicates that this fiux is most likely 
due to the primary star. We have convolved the lUE spectra with the GALEX FUV and 
NUV filter bandpasses in order to compute the GALEX-equivalent fiuxes (or upper limits) 
for these sources. We find, for R Lmi, 3cr upper limits of 0.032 mJy (FUV) and 0.16 mJy 
(NUV). For the other two sources, the derived fiuxes {3a errors) are as follows: TW Hor - 



•^taken from http://galexgi.gsfc.nasa.gov/docs/galex/tools/Resolution_Response/index.html 
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0.096 ± 0.051 mJy (FUV) and 0.41 ± 0.05 mJy (NUV); V Hya - 0.12 ±0.11 mJy (NUV). 

We now examine two plausible explanations for the FUV excesses, both of which involve 
the presence of a companion star. 



3.1. A Hot Companion 

The NUV and FUV excesses may result from the presence of a companion star which 
is significantly hotter than the primary. We have therefore made least-squares fits to the 
FUV and NUV excesses of each object by including the contribution of a companion star 
(Table[2]). We have used models by Kurucz (Castelli & Kurucz 2002) for the companion 
spectra. 

For RWBoo, AACam and VEri, the same value of Ay was applied to the model spec- 
trum of the companion, as derived from fitting the primary. In the case of R UMa, which is 
listed in GCVS with a spectral type M3-M9, the best fits to the optical and near-infrared 
fluxes were obtained with Ay = 2.6 (0.0) for an M3 (M9) primary spectrum. Therefore, in 
our least-squares fitting of the companion, we tried three values of Ay (0,1.3 & 2.6) to scale 
the companion black-body spectrum. 

Our modelling (Table[2]) provides the fractional luminosity (relative to the primary) of 
the companion. For RWBoo and VEri, the Hipparcos parallaxes, 3.09 ± 1.1 mas and 4.56 ± 
1.08 mas, give distances of 320 and 220 pc, respectively, implying companion luminosities of, 
Lc= 18 and 6 Lq for these two sources. For AA Cam and R UMa, the parallax measurements 
are not significant, and the luminosities given in Tableware for a nominal distance of 0.5 kpc. 

The value of Lc for RW Boo is consistent with that expected for a mid-A main-sequence 
star, but it is too low in the case of AA Cam, R UMa, and V Eri, since the luminosities 
of main-sequence stars using our most favorable (i.e. lowest) model values for each of 
these, Te// = 7800, 8900, 9250 K (i.e., spectral type ~A6 to ~Al-A2), lie in the range ~(10- 
35) Lq (Cox 2000, Table 15.7). An appeal to distance ambiguities for these 3 stars does not 
help to resolve the problem of the derived Lc values being too low for main-sequence stars. 
We have tried increasing the source distance in order to bring up L^ to its main-sequence 
values. For AACam and RUMa, this exercise results in Lp=2A x IO^Lq and 9.3 x 10'^Lq, 
using the most favorable models values in Table 2 - i.e., the lowest values of T^ff (7800 K 
and 8900 K) and the highest values of Lc (1.79 and 1.01 L©). The value of Lp for RUMa 
(AACam) is certainly (probably) too high for an AGB star. For VEri, a factor 2.4 increase 
in the distance is needed to scale up Lc to a main-sequence luminosity, but such a large in- 
crease is significantly beyond that allowed by its parallax data (factor 1.3). We rule out the 
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possibility that the companions are low- luminosity white dwarfs (WDs) on cooling tracks, 
because stellar evolutionary models show that by the time WDs have cooled to 10'* K, their 
luminosities are orders of magnitude below ILq. 

In our models we adopted the extinction curve as tabulated by Whittet (1992). Our 
quoted modeling uncertainties do not take into account uncertainties in the extinction curve 
at NUV and FUV wavelengths. We repeated the modeling using extinction curves for the 
LMC Supershell and the SMC Bar (Gordon et al. 2003), which along with the Galactic 
extinction roughly cover the range of curves found in circumstellar dust and are well-studied. 
The results are shown in Table [31 Although the best-fit temperatures all shifted up when the 
LMC and SMC extinction curves were used, the cooler companion models (i.e., for RWBoo 
and AA Cam) proved reasonably insensitive to the choice of extinction curve, but dramatic 
differences were found for RUMa and VEri. In general, we could not obtain a good fit 
even with the highest-temperature models available (39000 K) while using the LMC or SMC 
curves for these two sources. 



3.2. Accretion onto a Companion Star 

Five out of nine objects in Tabled] were observed on more than one epoch in one or both 
of the GALEX bands - in each instance, significant photometric variability was observed 
(Table [1]). We have checked that this variability is not due to systematic calibration uncer- 
tainties because the average and median fractional differences of the fluxes for the brightest 
40 field objects in the images from the different epochs are negligible. 

A plausible interpretation for the photometric variability is related to the presence of a 
nearby companion. This interpretation is motivated by ultraviolet observations of Mira, a 
symbiotic star in which Mira A is the AGE primary and MiraB is a compact companion (at 
a separation of &!6) which is accreting matter from MiraA's wind. lUE spectra of Mira in 
the wavelength region covered by the GALEX bands show the presence of strong emission 
lines ascribed to Mira B (Reimers & Cassatella 1985). The strongest of these (due to 0,N, & 
C, seen by lUE during 1979-95) were found to fade by a factor >20 by 1999-2001, and then 
start increasing back to their original levels by 2004 (Wood & Karovska 2006). Assuming 
the revised Hipparcos-based distance of 107 pc to Mira, the combined maximum fluxes of 
such emission lines, if present in our sources, would correspond to an artificial continuum in 
GALEX's broad-band FUV filter of about 0.1 mJy at a typical source distance of 350 pc, thus 
comparable to the measured values of the FUV fluxes in our sources. During the lUE era 
(1979-80 & 1990-95), the FUV lines and continuum varied by a factor ~2 in Mira. Karovska, 
Wood and co-authors (Wood, Karovska & Hack 2001, Wood, Karovska & Raymond 2002) 



conclude that the UV variabihty most hkely results from variations in the accretion rate 
onto MiraB. Although Mira's variability has been observed on a much longer time-scale 
than the ones sampled in our GALEX data, accretion of matter from the primary AGB 
wind onto a companion provides a plausible explanation for the presence of FUV emission 
and its variability in our sources. 



4. Discussion 

Our small survey of 21 AGB stars for UV excesses has resulted in a substantial number 
of NUV and/or FUV detections. Nine of these were detected in the FUV band and are the 
subject of this paper. A detectable FUV flux at even a few //Jy level is several orders of 
magnitude too high to be explained by photospheric emission from the relatively cool primary 
stars in our sample, and hence is an "excess" which requires an alternative explanation - 
most likely the presence of a binary companion. The excesses arise either as a result of 
photospheric emission from a hotter companion, and/or from an accretion disk around the 
companion. Spectroscopic monitoring in the FUV of these sources is needed in order to 
distinguish between these two mechanisms. 

We detected NUV fluxes in 19/21 of our objects with high S/N, many of which are 
also likely to be "excesses", but for which such an inference is more uncertain because of 
the significantly larger contribution of the primary in the NUV compared to the FUV. A 
discussion of the detection statistics of, and the biases in, our full sample, is deferred to a 
forthcoming paper. 

Although our discovery of a UV-excess attributable to a different star than the primary 
AGB star does not directly imply that the former is a gravitationally-bound companion, it 
is the most likely explanation. This is because the UV sky is rather "empty" (i.e., much 
more scarcely populated than at optical wavelengths), hence the probability, p false that the 
FUV-emitting object is simply positionally coincident on the sky with the primary (i.e., 
lying within a radius of 2" from the primary), is very small. Using the object number count 
(per deg^ mag) versus magnitude plot for the GALEX FUV band (Bianchi et al. 2007), 
we find that for AACam (the faintest of our modelled sources), p false for an object of FUV 
magnitude lying within a 0.5 mag bin centered around the FUV mag of AA Cam, is 8 x 10^^. 
The Galactic latitudes of our objects are similar to those of the fields used by Bianchi et al., 
hence it is appropriate to use their point-source densities. 

V Hya, which has the largest FUV fiux, as well as the highest FUV-to-NUV fiux ratio 
amongst all our targets, is well known for its collimated, high-velocity, outfiows, an extended 
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dusty torus, and an inner hot disk. The outflows were first seen via infrared absorption hues 
in the CO 4.6 ;um vibration-rotation band (Sahai & Wannier 1988); recent interferometric 
mapping of the miUimeter-wave CO hne emission shows the collimated structure of the fast 
outflow (e.g., Hirano et al. 2004). More recently, observations with STIS/HST reveal the 
presence of a high-velocity blob moving away from the central source at (projected) speeds 
upto 220km s^^, and a hot, slowly expanding (10-15 km s~^) central disk-like structure 
(Sahai et al. 2003). Although the expansive kinematics of the latter implies that it is not 
an accretion disk, it may result from a recent phase of equatorially enhanced mass-loss, 
which may be enhancing the accretion process. VHya is thus the best example to date of 
an evolved star with an active, collimated outflow, dense equatorially-flattened structures 
possibly related to a central accretion disk, and an inferred binary companion from our UV 
excess measurements. 

We would like to thank an anonymous referee for his/her thoughtful review of our 
paper. We acknowledge discussions with Patrick Morrissey related to the possibility of red 
leaks in the GALEX FUV and NUV bands. RS's contribution to the research described 
in this publication was carried out at the Jet Propulsion Laboratory, California Institute 
of Technology, under a contract with NASA. RS thanks NASA for financial support via 
a GALEX award and an LTSA award. KF was partially funded by a SURF scholarship 
and through the Cornell Presidential Research Scholars (GPRS) program.AGdP is partially 
financed by the Spanish Ramon y Cajal program and the Programa Nacional de Astronomia 
y Astrofisica under grant AYA 2006-02358. CSC is partially funded for this work by the 
Spanish MCyT under project AYA 2006-14876 and the Spanish MEC under project PIE 
2007501028. 
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Table 1: 


AGB Stars with UV Excesses 


Target 


Band 


Epoch'' 


Exp. Time (s) 


Flux (mjy) 


RW Boo'' 


FUV 


4220.75 


1726 


0.026 




NUV 


3861.15 


3396 


0.47 






4220.75 


1726 


0.39 






Average 


5122 


0.44 


AACam 


FUV 


3377.48 


1693 


0.014 




NUV 


3377.48 


1532 


0.22 






3425.29 


1693 


0.28 






Average 


3225 


0.25 


TDra 


FUV 


3984.2 


1580 


0.0055 




NUV 


3541.9 


1250 


0.017 






3587.5 


1309 


0.018 






3984.2 


1580 


0.032 






Average 


4151 


0.021 


VEri 


FUV 


3678.20 


1620 


0.060 




NUV 


3678.20 


1704 


0.14 


TWHor'' 


FUV 


3349.3 


4381 


0.026 






3351.2 


1345 


0.032 






3671.8 


109 


0.034 






3706.6 


255 


0.023 






3706.7 


586 


0.033 






4018.8 


747 


0.034 






Average 


7423 


0.027 




NUV 


3349.3 


4381 


0.53 






3351.2 


1345 


0.55 






3359.25 


2214 


0.34 






3671.8 


109 


0.78 






3706.6 


339 


0.26 






3706.7 


645 


0.27 






4018.8 


747 


0.99 






Average 


9780 


0.49 


VHya 


FUV 


3421.86 


1705 


0.12 






3778.54 


1003 


0.15 






Average 


2708 


0.13 




NUV 


3421.86 


1705 


0.11 






3778.54 


1003 


0.13 






Average 


2708 


0.12 


AF Peg'^ 


FUV 


3280.70 


1551 


0.011 




NUV 


3280.70 


1551 


0.090 


RUMa 


FUV 


3742.14 


1703 


0.041 




NUV 


3742.14 


1703 


0.12 


VYUMa 


FUV 


3742.50 


1704 


0.0061 




NUV 


3742.50 


1704 


0.23 



"JD-2450000 

''Observed at multiple epochs; data from epochs separated by <0.1 day are averaged together in the individual- 
epoch rows 

^Observed at 2 epochs in the NUV band only; quoted fluxes are for AFPeg using a 2-Gaussian fit to the 
image of the latter and the partially blended nearby star, for Epoch 1 
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Table 2: Model Resul 


ts 










Target 


Primary 


D 


Av 


TefS^ 


Lc"" 


-l^p/ -l^c 




Sp.Type 


(kpc) 




(K) 


(Le) 




RWBoo 


M5 


0.32 


2.3 


8200 (-500,300) 


18 (-5,7) 


280 (-80,110) 


AACam 


M5 


0.5 


1.0 


8200 (-400,400) 


1.1 (-0.6,0.7) 


3200 (-1300,4500) 


VEri 


M6 


0.22 


2.9 


10000(-700,4100) 


6.2 (-2.9,2.2) 


910 (-240,810) 


RUMa 


M3-M9 


0.5 


1.3 


9200(-300,+1100) 


0.85(+0.2,-0.4) 


5300(-900,+4700) 



"The numbers in parenthesis represent 3cr modclhng uncertainties; temperature and luminosity uncertainties 
are inversely correlated 



Table 3: Best-fit temperatures (in K) under different extinction curves. 



Target 


Av 


Galactic*^ 


LMC^ 


SMC^ 


RWBoo 


2.3 


8200 (-500,300) 


8700 (-500,400) 


No fit 


AACam 


1.0 


8200 (-400,400) 


8500 (-400,500) 


9000 (-400,1400) 


VEri 


2.9 


10000 (-750,4500) 


No fit 


No fit 


RUMa 


1.3 


9200(-300,+1100) 


33000 (-22000)'^ 


No fit 



"The numbers in parenthesis represent 3a modelling uncertainties. 

''The 3(7 upper limit for RUMa fell beyond the range of available models. 



- 12- 

REFERENCES 

Bergeat, J., Knapik, A., & Rutily, B. 2001, A&A, 369, 178 
Bianchi, L., et al. 2007, ApJS, 173, 659 

Bodenheimer, P., Burkert, A., Klein, R. I., & Boss, A. P. 2000, Protostars and Planets IV, 

675 

Boggess, A., et al. 1978a, Nature, 275, 372 

Boggess, A., et al. 1978b, Nature, 275, 377 

Bond, H.E. 2000, ASP Conf. Ser., Vol. 199. Eds. J.H. Kastner, N. Soker, & S. Rappaport p. 
115 

Castelli, F. & Kurucz, R.L. 2003, in "Modelling of Stellar Atmospheres", lAU Symposium 
210, eds. by Piskunov, W.W. Weiss, and D.F. Gray, ASP Conf. Ser., p. A20 

Ciardullo, R., M. S., Fullton, L. K., Zhang, C.-Y., & Schaefer, K. G. 1999, AJ, 118, 488 

Corradi, R.L.M., & Schwarz, H.E. 1995, A&A, 293, 871 

Cox, A.N. 2000, Allen's Astrophysical Quantities, 4th ed., AIP Press: Springer 

De Marco, O., Bond, H.E., Harmer, D., & Fleming, A.J. 2004, ApJ, 602, L93 

Duquennoy, A. & Mayor, M. 1991, A&A, 248, 485 

Garcia-Segura, G. 1997, ApJ, 489, L189 

Fluks, M. A., Plez, B., The, P. S., de Winter, D., Westerlund, B. E., & Steenman, H. C. 
1994, A&AS, 105, 311 

Gordon, K. D., Clayton, G. C, Misselt, K. A., Landolt, A. U., & Wolff, M. J. 2003, ApJ, 
594, 279 

Hirano, N., et al. 2004, ApJ, 616, L43 

Knapp, G. R., Pourbaix, D., Platais, I., & Jorissen, A. 2003, A&A, 403, 993 

Rich, R. M. 2005, ESO astrophysics symposia: "Multiwavelength Mapping of Galaxy For- 
mation and Evolution", Eds. A. Renzini & R. Bender, (Springer: Berlin), 203 

Riemers, D. & Cassatella, A. 1985, ApJ, 297, 275 



-13- 

Samus, N. N., Durlevich, O. V., & et al. 2004, VizieR Online Data Catalog, 2250, 

Sahai, R. k Trauger, J.T. 1998, AJ, 116, 1357 

Sahai, R., & Wannier, R G. 1988, A&A, 201, L9 

Schwarz, H.E., Corradi, R.L.M. & Melnick, J. 1992, A&A Suppl.Ser., 96, 23 

Soker, N. 1998, ApJ, 496, 833 

Soker, N. & Rappaport, S. 2001, ApJ, 557, 256 

Sorensen, P. M., & Pollacco, D. L. 2003, Astronomical Society of the Pacific Conference 
Series, 303, 494 

van Winckel et al. 1999, A&A, 343, 202 



, ^ ilJO-l i, ^ J.^, 



Whittet, D. C. B. 1992, Dust in the Galactic Environment, Institute of Physics Publishing, 
p.67 

Wood, B. E., Karovska, M., & Hack, W. 2001, ApJ, 556, L51 

Wood, B. E., Karovska, M., & Raymond, J. C. 2002, ApJ, 575, 1057 

Wood, B. E., & Karovska, M. 2004, ApJ, 601, 502 

Wood, B. E., & Karovska, M. 2006, ApJ, 649, 410 



This preprint was prepared with the AAS I^TJtjX macros v5.2. 



-14- 



o 



s I 



X 
3 

o 



I 



-I— I— I— p-i— I— I— I— [■ 
RWBoo 




o 



Ln 
I 



o - 



Ln 
I 




0.5 1 1.5 

Wavelength (micron) 



0.5 1 1.5 2 

Wavelength (micron) 



Fig. 1. — GALEX (NUV & FUV) and ground-based (optical & near-IR) fluxes [red sym- 
bols) of AA Cam, V Eri & RUMa, with model spectra {blue: cool AGB star, green: hot 
companion). The expected NUV flux {blue/ cyan diamonds) from the cool AGB star (the 
expected FUV flux due to AGB star lies below the minimum of the flux range), and from 
the hot companion + cool AGB star [green diamonds) are also shown. The inset shows an 
expanded view of the UV-Blue region. The model NUV/FUV fluxes {diamond symbols) are 
obtained by convolving the model spectrum with the GALEX NUV and FUV bandpasses 
{black curves, inset). 



